NUWC-NPT  Technical  Report  10,817 
19  January  1995 


Preface 


This  work  was  conducted  under  the  Submarine  Communications  Exploratory  Development 
Project  (SUBCOMMS),  RC32C16.  This  work  was  conducted  in  support  of  the  "Low  Profile 
Antenna"  task  under  the  SUBCOMMS  Project.  The  Principal  Investigator  for  this  task  is  Keith 
Lewis,  Code  3433.  The  SUBCOMMS  Project  manager  is  Frederick  Allard,  Code  3496.  TTie 
SUBCOMMS  Project  sponsor  is  Dr.  Sherman  Gee,  ONR  313. 

The  technical  reviewer  for  this  report  was  Kurt  F.  Hafner  (Code  3413). 


Acknowledgments 

The  authors  wish  to  thank  Mr.  Kurt  Hafner  of  the  Communications  Antennas  Branch  and  Dr. 
Rajeev  Bansal  of  the  University  of  Connecticut  for  their  careful  review  of  this  report.  In  addition, 
the  authors  wish  to  thank  Mr.  Mark  Casati  of  the  Antisubmarine  Warfare  Systems  D^artment  for 
his  assistance  in  obtaining  the  computed  results. 


Reviewed  and  Approved:  19  January  1995 


D.  M.  Viccione 

Head,  Submarine  Electromagnetic  Systems  Department 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  fopofting  burden  kx  mis  calectioo  c<  informaiion  «  wtimaed  to  awefage  1  hour  p«r  resoonsa,  including  the  tinw  for  rw»v»^  iminxaton*,  s«»c*iing  exiting  data  soucas.  ^attwing  and 
maintainjnfl  tt>«  data  needed,  arid  completing  and  reviewing  the  collection  of  information.  Sand  comments  regarding  this  burden  estimate  or  other  aspect  of  the  collection  of  infomeiion. 
Including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  lor  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway  Suite  1204* 
on.  VA  22202-4TO  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704^1881.  W»hinqton.  DC  20503. 


3.  REPORT  TYPE  AND  DATES  COVERED 


Final 


5.  FUNDING  NUMBERS 


1.  AGENCY  use  ONLY  (Leave  Blank) 

Z  REPORT  DATE 

19  January  95 

4.  TITLE  AND  SUBTITLE 

Approximate  Capacitance  Formulas  for  Electrically 
Small  Tubular  Monopole  Antennas 


6.  AUTHOR{S) 

David  F,  Rivera,  Dr.  John  P.  Casey 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Undersea  Warfare  Center  Detachment 
39  Smith  Street 

New  London,  Connecticut  06320-5594 


a  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TR- 10,817 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 
Code  313 
Washington,  DC 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


1 2a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  wofds) 

Analytical  approximations  have  been  developed  to  compute  the  capacitance 
of  electrically  small  tubular  monopole  antennas.  The  expressions  developed 
are  algebraically  simple,  facilitating  manual  calculations  and  rapid  iteration 
of  designs  when  used  in  a  computer  program.  Plots  are  given  of  the  accuracy 
of  previous  expressions  and  the  approximations  developed  herein  as  a  function 
of  the  arguments.  It  is  shown  that  the  approximations  are  sufficiently 
accurate,  making  them  useful  as  tools  to  aid  in  the  design  of  electrically 
short  cylindrical  monopoles  over  a  wide  range  of  heights,  lengths,  and 
diameters . 


14.  SUBJECT  TERMS 

Electrically  small  antennas 
Monopoles 


Capacitance 
Moment  Method 


15.  NUMBER  OF  PAGES 
70 


16.  PRICE  CODE 


17:  SECURITY 

CLASSIFICATION  OF 
REPORT 

UNCLASSIFIED 


ia  SECURITY 
.  CLASSIFICATION  OF  THIS 
PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION  20.  UMITATION  OF  ABSTRACT 
OF  ABSTRACT 


UNCLASSIFIED 


UNLIMITED 


Standard  Form  298  (Rev  2-89) 

Prascribed  by  ANSI  Std  Z^3Q■^^ 


TR  10,817 


TABLE  OF  CONTENTS 

Section  Page 

LIST  OF  ILLUSTRATIONS .  u 

LIST  OF  TABLES .  iv 

1.  BACKGROUND . 1 

2.  APPROACH .  2 

3 .  EXISTING  APPROXIMATIONS  AND  THEIR  REGIONS  OF  VALIDITY ...  3 

3.1  Thin  Tubes .  3 

3.2  Thick  Tubes .  5 

4 .  EXTENSION  OF  GROVER’S  FORMULA .  6 

5.  APPROXIMATE  EXPRESSION  FOR  INTERMEDIATE  PARAMETER 

RANGE . 8 

6 .  DISCUSSION  OF  RESULTS . 10 

7.  APPUCATION  TO  A  TUBE  OF  NONCmCULAR  CROSS  SECTION .  11 

8.  SUMMARY.. . 11 

APraNDIX  A  -  CAPAQTANCE  FORMULAS  FOR  TUBULAR 

MONOPOLE  FEEDS  AND  OTHER  EFFECTS .  25 

1.  Introduction  .  25 

2.  Summary  of  Capacitance  Formulas . 25 

APPENDIX  B  -  APPROXIMATION  FOR  THE  RATIO  OF  THE 

COMPLETE  ELLIPTIC  INTEGRALS  KlkO/Kfk)  .  45 

1.  Background .  45 

2.  Derivation  of  K(kO/K(k)  Approximation  . 46 

3.  Discussion  of  Results . 49 

APPENDIX  C  -  NUMERICAL  CONVERGENCE  OF  THE  METHOD 

OF  MOMENTS  ALGORITHM .  57 

1.  Introduction . 57 

2.  Results .  57 

REFERENCES.... .  65 


1 


TR  10,817 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  (a)  Tubular  monopole. 

(b)  Coordinate  system  used  for  tubular  monopole  analysis .  12 

2  Typical  methods  of  feeding  tubular  monopoles. 

(a)  Wire  feed,  (b)  Conical  feed .  13 

3  Charge  distribution  along  a  tubular  monopole .  14 

4  Error  contours  for  Grover's  formula  (6) .  15 

5  Model  used  in  the  derivation  of  the  conformal  mapping  approximation  (8). 

(a)  Tubular  dipole,  (b)  Equivalent  coplanar  stripline .  16 

6  Comparison  of  normalized  capacitance  using  CMA  formula  (8)  and 
measured  data.  The  dashed  lines  are  the  error  bounds  in  the 

measurements . 17 

7  Error  contours  for  the  conformal  mapping  approximation  (8) .  18 

8  Error  contours  for  the  extended  Grover  formula  (17) .  19 

9  Surfece  of  the  normalized  tubular  capacitance  used  to  determine  ACF  (23).....  20 

10  Errw  contours  for  the  approximate  capacitance  formula  (23) .  21 

1 1  Regions  of  validity  (oror  <  10% )  of  tubular  monopole  capacitance  formulas . .  22 

A- 1  Semi-infinite  coplanar  plate  cross  section. 

(a)  Symmetrical  arrangemeit;  (b)  unsymmetrical  equivalent  .  32 

A-2  Thickness  correction  for  symmetrically  arranged  coplanar  plates .  33 

A-3  Laigthwise  view  of  a  capacitive  iris  in  a  waveguide.  The  nearest 
extraneous  surface  in  the  illustration  is  the  waveguide  wall. 

From  Cohn  [26] . 34 

A-4  Thin-wire  feed  configuration .  35 

A-5  Cross  section  of  coaxial  line/monopole  transition  (w=0).  The 

quantities  2a  and  2b  denote  the  inner  and  outer  diameters,  respectively .  36 

A-6  The  normalized  transition  capacitance  for  a  monopole  driven  by  a 

coaxial  line.  Figure  from  King  [27] . . .  37 

A-7  Truncated  cone  with  tip  pointed  toward  the  ground  plane. 

Tip  is  close  to,  but  does  not  touch  the  ground  plane .  38 

A-8  Plate  of  arbitrary  shape  over  a  ground  plane  with  interposed  dielectric. 

Thickness  (or  height)  of  dielectric  is  h  .  39 


11 


TR  10,817 


LIST  OF  ILLUSTRATIONS  (Cont’d) 

Figure  Page 

A-9  Arrangement  of  variables  used  in  computing  the  function  Ai(P) 

in  (A- 1 3a).  From  Kuester  [18]  . . . . . .  40 

B- 1  Ratio  of  the  complete  elliptic  integrals  of  the  first  kind,  K(k0/K(k); 

solid  line.  Shown  in  dotted  lines  are  the  limiting  values  of  K(k')/K(k), 
obtained  from  (B-3) .  50 

B-2  Cross  section  of  a  curved  plate  transmission  line  used  to  derive  the 
first  approximation,  (B-6).  The  angular  separation,  0,  is  related  to 
s  and  D  by  sin  9  =  s/D . 51 

B-3  Values  of  5  required  to  yield  correct  values  of  KQcO/KOc) .  52 

B-4  Relative  error  incurred  in  computing  KOcO/KOc),  using  (B-7) .  53 


B-5  Value  of  (B-7)  relative  to  (B-3b)  in  the  region  of  poorest  accuracy  (k  1) . . . .  54 


C- 1  Normalized  capacitance  of  a  tubular  monopole  versus  the  number  of 
basis  functions  for  logio(D)  =  0  and  logio(H)  =  -1  as  determined  by 
the  method  of  moments  algorithm .  59 

C-2  Normalized  capacitance  of  a  tubular  monopole  versus  the  number  of 
basis  functions  for  logio(D)  =  0  and  logio(H)  =  -3  as  determined  by 
the  method  of  moments  algorithm.... .  60 

C-3  Number  of  basis  functions  required  to  insure  a  minimum  of  three 

significant  figure  accuracy  in  the  method  of  momaits  algorithm 
used  for  computing  tubular  monopole  capacitance .  61 

C-4  Variation  of  the  normalized  tubular  monopole  capacitance  with  height 

forlogio(D)  =  -l .  62 

C-5  V ariation  of  the  normalized  tubular  monopole  capacitance  with  height 

for  logio(D)  =  3 . . .  63 


TR  10,817 


LIST  OF  TABLES 

Table  Page 

I  Ranges  of  validity  of  tubular  capacitance  formulas  (error  <  10%) .  23 

A-1  Comparison  of  coaxial  end-effect  approximation  (A-9)  with  theory  and 

experiment . 41 

A-2  Values  of  the  coefficient  Ai(P)  used  with  capacitance  approximations 

(A- 12)  and  (A- 13b)  for  various  plate  geometries .  42 

A-3  Effective  radius  (ag)  used  with  capacitance  approximations  (A- 14) 

and  (A- 15)  for  various  plate  geometries .  43 

B-1  Comparison  of  the  approximation  (B-7)  with  exact  values,  for 

modulus  k  close  to  unity . 55 


IV 


TR  10,817 


1.  BACKGROUND 

A  monopole  antenna  is  said  to  be  electrically  small  when  its  largest  physical  dimension  is  much 
smaller  than  the  wavelength  at  which  it  operates.-  These  antennas  are  commonly  found  in  the 
portion  of  the  radio  ftequency  spectrum  spanning  from  VLF  to  MF  (3  -  3000  kHz).  At  these 
frequencies,  such  antennas  are  employed  in  maritime  communications,  air  and  coastal  navigation, 
as  well  as  local  broadcasting  [  1  -3] . 

The  input  impedance  of  an  electrically  small  monopole  can  be  represented  by  a  series  circuit 
comprising  a  resistance  R  and  a  capacitive  reactance  -l/coC,  such  that  R  «  1/coC  [1].  The  series 
resistance  is  the  sum  of  the  radiation  resistance  and  other  resistances  attributed  to  ohmic 
(dissipative)  losses.  The  monopole  capacitance  C  can  be  derived  using  electrostatic  methods  and 
depends  on  three  principal  quantities:  its  diameter  d,  length  I,  and  height  h  measured  from  the 
bottom  of  the  antenna  to  the  ground  plane.  Such  an  arrangement  is  shown  in  Fig.  1(a),  where  the 
wall  thickness  t  of  the  tube  is  assumed  to  be  very  thin  compared  to  the  diameter. 

In  the  design  of  such  an  antenna,  the  goal  is  to  minimize  the  capacitive  reactance  and  maximize 
the  radiation  resistance.  Achieving  this  goal  usually  involves  an  increase  in  the  effective  antenna 
radius  and/or  length  by  using,  for  example,  a  folded  monopole  with  one  or  more  input  impedance 
transforming  elements  [4-5],  a  capacitive  top  load  [6],  or  both  [7].  In  addition,  the  capacitance  of 
the  electrically  small  tubular  monopole  is  extremely  important  for  the  determination  of  the  antenna's 
power  handling  ability  (P^ax)  and  the  bandwidth-radiation  efficiency  product  (BW  t|)  [8-9]  given 
as  follows: 


P 


max 


_640K^f*Y^hiC'^ 

r2 


(1) 


BW  Ti 


320  7l3f^he^C 
c2 


(2) 


where 

hg  =  effective  height, 
f  =  frequency, 

C  =  antenna  capacitance, 

Vj,  =  maximum  allowable  base  voltage  before  breakdown, 
c  =  velocity  of  light. 

In  designing  electrically  small  monopole  antennas,  the  engineer  must  be  concerned  with  these 
antenna  performance  parameters  as  well  as  any  number  of  mechanical  requirements  that  may  be 
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imposed.  The  complete  design  process  is  therefore  seen  as  highly  iterative,  involving  the 
examination  of  achievable  electrical  and  mechanical  properties  until  a  satisfactory  solution  is 
obtained. 

A  precise  determination  of  the  input  capacitance  of  an  electrically  small  antenna  may  be 
obtained  through  the  solution  of  the  potential  integral  equation  for  the  unknown  charge  distribution 
using  the  method  of  moments  [10].  However,  such  a  numerical  method  does  not  lend  itself  to 
rapid  itaative  design  calculations. 

Approximate  analytical  expressions  for  the  input  capacitance  are  available  for  two  extreme 
cases,  i.e.,  for  a  thin  wire  (d//  «  1)  and  for  a  thick  tube  [(d// »  1),  with  vanishingly  thin  walls], 
each  above  a  ground  plane.  The  capacitance  of  a  thin  vertical  wire  above  a  ground  plane  was 
d^ermined  by  Grover  [1 1].  By  the  use  of  a  method  originally  proposed  by  Howe  [12],  GrovCT 
obtained  an  analytical  expression  for  the  capacitance  by  assuming  that  the  charge  distribution  along 
the  antenna  was  constant.  The  converse  problem,  that  of  the  capacitance  of  a  thick  tubular 
monopole,  was  solved  by  Casey  and  Bansal  [13].  Through  modification  of  the  per-unit  loigth 
capacitance  of  a  coplanar  stripline  given  by  Hanna  [14],  Cas^  and  Bansal  obtained  an  expression 
for  the  equivaloit  tubular  monopole  in  terms  of  elliptic  integrals. 

Although  the  expressions  died  above  have  been  experimaitally  validated  for  some  spedfic 
monopole  dimensions,  the  extait  of  diameters,  lengths  and  ground  plane  separations  for  which 
they  are  accurate  has  been  unknown.  When  used  for  the  rapid  iteration  of  a  design,  capacitance 
values  obtained  with  the  known  expressions  will  be  of  questionable  value.  A  determination  of  the 
regions  of  validity  of  the  existing  tubular  monopole  capacitance  formulas  is  thwefore  necessary  for 
their  successful  implementation.  It  is  the  purpose  of  this  rqx)rt  to  summarize  the  results  of  such  an 
investigation  and  present  new  capadtance  formulas  for  param^a*  ranges  not  covered  by  the 
existing  expressions.  The  end  result  of  this  study  is  a  collection  of  formulas  that,  taken  together, 
extend  the  parameter  range  so  that  the  computation  of  capacitance  can  be  facilitated  with  reasonable 
accuracy  for  almost  all  design  situations  considered  in  practice. 

2.  APPROACH 

The  range  of  validity  of  Grover’s  formula  for  the  capadtance  of  a  thin  wire  or  tubular 
monopole  and  the  range  of  validity  of  the  thick  tubular  monopole  formula  developed  by  Casey  and 
Bansal  will  be  examined  by  comparison  with  results  obtained  by  a  method  of  moments  solution 
[15].  (It  may  be  noted  that  the  method  of  momaits  results  have  been  experimentally  verified 
[13,15].)  The  range  of  validity  is  defined  here  as  that  area  where  the  formula  agrees  to  within  10% 
of  the  method  of  moments  results.  A  10%  error  region  is  chosen  since  influences  such  as  the 
antenna's  proximity  to  otha:  objects  or  irregularities  in  the  ground  plane  may  introduce  variations 
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of  this  magnitude  in  the  observed  capacitance. 

With  the  regions  of  validity  defined  for  the  existing  capacitance  formulas,  additional 
expressions  will  be  presented  that  are  suitable  for  parameter  ranges  not  covCTed  by  the  existing 
equations.  An  error  analysis  of  these  new  capacitance  formulas  will  follow. 

Tubular  monopoles  can  be  fed  in  a  variety  of  ways,  two  of  which  are  shown  in  Fig.  2.  In 
genml,  the  monopole  capacitance  is  approximately  the  sum  of  the  individual  capacitances  of  both 
the  tube  and  the  feed  sections.  The  capacitance  of  a  thin  feed  wire  (Fig.  2(a))  is  normally  much 
smaller  than  the  tube  capacitance  and  can  be  neglected.  In  contrast,  the  capacitance  of  the  conical 
feed  as  shown  in  Fig.  2(b)  must  usually  be  accounted  for.  Information  on  the  computation  of  the 
illustrated  feed  capacitances  can  be  found  in  [16- 17].  In  the  case  where  there  is  a  plate  attached  to 
both  the  feed  wire  and  the  lower  aid  of  the  tube,  a  first  order  estimate  of  the  monopole  capacitance 
is  the  sum  of  the  individual  capacitances.  Information  on  the  capacitance  of  a  flat  plate  is  given  in 
[18].  The  effects  of  top  loading  and  base  supports  may  be  treated  by  techniques  given  by  Bekose 
[3]  and  are  not  consido'ed  hare.  Appendix  A  provides  a  listing  of  capacitance  formulas  for  various 
tubular  monopole  feeds  along  with  otho*  effects  including  the  finite  wall  thickness  of  the  tube 


3.  EXISTING  APPROXIMATIONS  AND  THEIR  REGIONS  OF  VALIDITY 


3.1.  Thin  Tubes 

Consider  the  tubular  monopole  and  associated  coordinate  system  as  described  in  Fig.  1(b). 

The  electrostatic  potential  at  any  point  along  the  tube  surface  due  to  an  axisymmetric  surface  charge 
doisity  a(z)  =  q(z)/7id  induced  on  the  tube  is  found  by  summing  the  contribution  of  the  charge 
along  the  cylinder  as  [15] 


where 


zO  -  K(z  +  zO]  dz'  , 


z  e  (h,h+0 


/•" 


(3) 


(4) 


In  (3),  q(z)  is  the  charge  per  unit  length  while  K(z  -  z')  and  K(z  +  z')  are  the  kernels  associated 
with  the  tube  and  its  image,  respectively.  Note  that  a  cylindrical  coordinate  system  is  used  in  (3) 
and  (4),  where  (p  doiotes  the  azimuthal  variable  while  the  primed  and  unprimed  coordinates  refo"  to 
the  source  and  observation  points,  respectively.  Since  the  tube  is  highly  conducting,  it  will  be  an 
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equipotential  surface  with  V(z)  =  V  (a  constant),  where  V  is  the  potential  of  the  tube  with  respect  to 
ground. 

The  charge  distribution  that  results  from  solving  the  integral  equation  (3)  when  the  potential 
V(z)  is  constant  was  calculated  using  the  method  of  moments.  A  representative  charge  distribution 
is  shown  in  Fig.  3.  The  charge  density  is  fairly  constant  across  the  interior  of  the  interval  (h,h+/) 
with  singularities  at  the  edges.  The  charge  density  is  asymmetric  about  its  midpoint  as  more  charge 
accumulates  on  the  ground  plane  side  of  the  tube.  This  asymmetry  becomes  more  apparait  as  the 
separation  h  decreases. 

For  thin  tubes  or  wires  (d// «  1),  Grover  derived  a  formula  for  the  capacitance  of  a  monopole 
using  Howe's  m^hod  of  approximation.  In  Howe’s  approximation  (contrary  to  the  physical 
reality)  the  total  charge  Q  is  assumed  to  be  uniformly  distributed  ovct  the  length  of  the  tube, 
thereby  reducing  (3)  to 


-zO-KCz  +  zOjdz'  , 


Z6  (h,h+/)  . 


(5) 


Note  that  the  potential  becomes  a  function  of  position  along  the  tube.  The  capacitance  is  estimated 
as  C  =  QA^av*  whCTe  V^v  is  the  average  of  V(z)  over  (h,h+/).  Upon  application  of  the  above 
procedure  followed  by  furthCT  simplification  under  the  condition  d//  «  1,  GrovCT  derived  the 
following  [11]: 


whCTe  D  =  d//  and  y  is  defined  as 


(6) 


Y=l  +  (l+H)ln(l+H)-(l+2H)ln(l  +  2H)  +  Hln(4H)  (7) 

where  H  =  h//.  For  simplicity,  we  will  adopt  the  symbols  D  and  H  throughout  the  remainder  of  the 
report.  Note  that  as  H  7  ->  1  -  In  2  and  Grover's  formula  yields  the  capacitance  of  a  tube  of 
length  I  in  free  space.  However,  as  H  0, 7  1  and  the  capacitance  does  not  diverge  as 

expected,  but  instead  approaches  a  constant. 

A  contour  plot  illustrating  the  percentage  error  of  the  capacitance  per  unit  laigth  of  Grover's 
formula  based  on  a  comparison  with  a  method  of  moments  solution  is  given  in  Fig.  4  for  -4  < 
iogio(H)  -  1  ’3  ^  logio(D)  <  0.  GrovCT’s  formula  produces  data  within  10%  of  the  moment- 

method  results  for  D  <  0.008  over  the  entire  range  of  ground  plane  sqparations  considared  in  this 
stutfy.  Although  Grover’s  formula  does  not  produce  the  correct  result  in  the  limiting  case  as  H 
0,  it  still  yields  sufficiently  accurate  data  for  the  small  values  of  H  plotted  with  D  <  0.008.  As  H 


4 


TR  10,817 


increases,  Grover's  formula  generally  improves  until  H  =  0.04,  where  the  contours  begin  to  level 
off.  Grover's  formula,  though  not  based  on  assumptions  that  accurately  describe  the  charge  and 
voltage  distributions,  provides  a  useful  expression  for  the  capacitance  within  the  regions  stated 
above. 

The  method  of  moments  code  that  was  used  in  this  comparison  is  based  on  the  solution  of  the 
potoitial  integral  equation  (3)  for  the  unknown  charge  density  q(z).  The  numerical  procedure 
employs  pulse  expansion  functions  and  point  matching.  The  kernel  K(0  of  the  potential  integral 
equation,  defined  in  (4),  is  expressible  in  terms  of  the  complete  elliptic  integral  of  the  first  kind,  for 
which  polynomial  approximations  exist  [19].  For  each  data  point  computed,  a  sufficient  numbo" 
of  basis  functions  was  chosen  to  ensure  convCTgaice  to  at  least  three  significant  figures.  The 
details  of  this  program  are  provided  in  [15].  An  examination  of  the  numerical  convOTgence  of  the 
method  of  momaits  algorithm  used  in  this  report  is  given  in  Appendix  C. 


3.2.  Thick  Tubes 

Casey  and  Bansal  [13]  developed  an  expression  for  the  capacitance  of  a  tubular  monopole 
through  comparison  with  a  formula  for  the  capacitance  of  a  coplanar  stripline.  Figure  5  illustrates 
the  approximate  equivalence  utilized  in  this  development.  The  expression  is  based  on  a  conformal 
mapping  and  is  given  by 


where 


and 


C  =  2jteod 


K(kO 

K(k) 


(8) 


k  =  -ii- 

^  1+H 


k'=fr^ . 


K(k)  is  the  compile  dliptic  integral  of  the  first  kind.  An  accurate  approximation  for  the  ratio  of  the 
complete  elliptic  integrals  in  (8)  is  givai  as 


K(k)  -  I  [  k  4(1.^-) 


(11) 


The  development  and  error  analysis  of  (1 1)  are  given  in  Appendix  B.  Equation  (8)  will  be  referred 
to  as  the  conformal  mapping  approximation  (CMA)  formula.  Figure  6  illustrates  a  comparison 
betweai  the  CMA  formula  and  experimental  data  for  tubes  of  various  dimensions  and  ground  plane 
separations.  It  can  be  seen  that  the  data  are  in  good  agreement. 
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A  contour  plot  illustrating  the  percentage  error  of  the  capacitance  per  unit  length  of  the  CMA 
formula  based  on  a  comparison  with  the  method  of  moments  solution  is  given  in  Fig.  7  for  -4  < 
logio(H)  <  1  and  -2  <  logio(D)  <  3.  The  CMA  formula  produces  data  within  10%  of  the  method 
of  moments  results  for  thick  monopoles  within  the  region  D  >  2/ln[l  +  (3/H)]  and  H  >  lO"'*.  The 

former  boundary  was  obtained  by  fitting  a  curve  to  the  calculated  boundary  points.  For 
sufficiently  thick  monopoles  (i.e.,  D  >  10),  the  error  contours  are  nearly  vertical  lines  (independent 
of  D)  since  the  differences  in  the  electric  field  lines  associated  with  the  monopole  and  the 
corresponding  stripline  change  very  little  with  increasing  tube  diameter. 

For  D  >  10,  the  error  in  the  CMA  formula  slowly  oscillates  with  separation  from  the  ground 
plane  over  the  entire  range  of  H  shown,  and  then  monotonically  increases  for  large  H  (not  shown 
in  Fig.  7).  The  increasing  error  in  the  (2MA  formula  for  large  H  occurs  because  the  electric  field 
lines  extending  from  the  inside  surface  of  the  tube  are  no  longer  the  same  shape  as  those  extending 
from  the  outside  surface.  The  required  symmetry  is  preserved  for  greater  ground-plane  separations 
when  the  diameter  of  the  tube  is  larger. 


4.  EXTENSION  OF  GROVER’S  FORMULA 


As  discussed  in  the  previous  section,  Grover's  capacitance  formula,  based  on  Howe's  method 
of  approximation,  is  valid  only  for  thin  tubes  while  the  CMA  formula  is  valid  for  thick  tubes.  In 
an  attempt  to  bridge  the  gap  between  the  CMA  and  Grover  formulas,  the  authors  have  extended  the 
Grover  formula  to  include  fewer  restrictions  on  D  and  H.  An  outline  of  the  development  of  this 
formula  along  with  results  are  given  below. 

Consider  again  expression  (3)  for  the  electrostatic  potential  along  the  surface  of  the  tubular 
monopole  due  to  a  surface  charge  density  q(z)/7id.  Since  the  surface  charge  distribution  is 
axisymmetric,  the  kernel  K(Q  defined  in  (4)  may  be  approximated  as 


(12) 


The  above  approximation,  sometimes  referred  to  as  the  reduced  kernel,  is  applicable  to  thin  tubes. 
Next  we  will  apply  Howe's  approximation  to  the  determination  of  the  tubular  monopole 
capacitance  using  the  reduced  kernel  (12).  The  substitution  of  (12)  into  the  potential  integral  (5) 
yields 
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V(z)  = 


47ceo/ 


(z  -  z')  +  (d/2)2  V  (z  +  z')^  +  (d/2)2 


dz  ,  z  e  (h,  h+0  (13) 


where  Q  is  the  total  charge  along  the  monopole.  The  evaluation  of  the  integral  in  (13)  yields 

W(7)  =  Q  /in  (z-h)  +  V(z-h)2H-(d/2)^  1  _  Uz  +  h  +  I)  +  ^/(z  +  h  +  l)^  +  (d/2)^l| 

47ieo/  \  [(z-h-O  +V(z-h-/)^  +  (d/2)2j  [_  (z  +  h)  +  V(z  +  h)2  +  (d/2)2  J| 

Note  that  the  final  result  follows  since  sinh'^  (x)  =  In  (x  +  Vx^  +  1 ).  The  average  potential  along 
the  monopole  is 


Jrh+/ 


where 


'P(D,H)  =  sinh-i(^j  -  (1  +  H)  sinh'ij^^^^j  +  (1  +  2H)  sinh-^(^-^-^^j  -  H  sinh-i(^) 

+  a .  Vl  +  (D/2)2  +  Vh2  +  (D/4)2  +  V(1  +  H)2  +  (D/4)2  -  V(1  +  2H)2  +  (D/4)2  (16) 


with  D  =  d//  and  H  =  h/l.  The  monopole  capacitance  is  given  by 

C=;^=  (17) 

^(D,H)- 

The  extended  GrovCT  formula  (17)  is  considerably  more  involved  than  Grover’s  formula  (6).  For 
the  case  of  very  thin  tubes  (D  «  1)  it  can  be  shown  that  (17)  reduces  to  (6). 

A  contour  plot  illustrating  the  pacentage  error  of  the  capacitance  per  unit  length  of  the  extended 
Grover  formula  based  on  a  comparison  with  the  method  of  moments  data  is  given  in  Fig.  8  for  -4 
<  logio(H)  <  1  and  -3  <  logio(D)  ^  1.  A  comparison  of  Figs.  4  and  8  indicates  that  the  extended 
Grover  formula  provides  a  small  improvement  over  Grover's  formula  for  H  >  0.1.  More 
specifically,  for  H  >  0.1,  one  is  able  to  model  monopoles  (within  10%  of  the  moment-method 
data)  with  the  extended  Grover  formula  for  D  <  1.0  in  comparison  to  D  ^  0.35  with  Grover’s 
formula.  Both  Grover’s  formula  and  the  extended  Grover  formula  produce  similar  results  for  H  < 
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0.1. 

The  inability  of  the  extended  Grover  formula  to  model  much  thicker  tubes  is  attributed  to  the 
approximate  reduced  kamel  (12)  which  is  valid  only  for  thin  tubes.  In  addition,  the  Howe 
approximation  limits  the  ability  of  the  extended  Grover's  formula  to  yield  an  improvemait  over 
Grover’s  formula  for  small  relative  ground  plane  separations  (H  <  0.1).  In  conclusion,  with  the 
added  number  of  terms  involved,  the  extended  Grover  formula  does  not  appear  to  offer  a 
significant  improvement  ova"  Grover's  result. 

5.  APPROXIMATE  EXPRESSION  FOR  INTERMEDIATE  PARAMETER  RANGE 

In  the  previous  sections,  the  capacitance  formulas  presaited  were  for  cases  in  which  the 
antenna  is  considered  thin  or  thick.  As  a  result,  an  expression  was  sought  that  produced  a  fit  to  the 
capacitance  data  obtained  from  a  momait-method  calculation  for  regions  not  accurately  represented 
by  either  the  CMA  or  Grover’s  formula.  The  area  of  interest  is  a  rectangular  region  with 
boundaries  defined  by  -2.5  <  logio(D)  ^  1  and  -4  <  logio(H)  <  1.  This  region  was  chosen  since  it 
was  considered  to  cover  most  of  the  areas  where  the  existing  formulas  fail  and  provides  a  sufficient 
amount  of  ovalap.  A  suitable  form  for  an  expression  that  adequately  describes  the  capacitance 
variation  with  variables  D  and  H  is  arrived  at  by  the  observations  that  follow. 

To  be  useful,  the  desired  capacitance  expression  must  possess  the  limiting  behaviors  described 
below. 


i)  Tube  close  to  the  groundplane 

In  this  region,  the  capacitance  can  be  represented  by  the  behavior  of  the  CMA  (8)  for  H  0. 
Here  one  can  utilize  the  asymptotic  rq)resentation  for  the  ratio  of  elliptic  integrals  K(kO/K(k)  for 
small  modulus  k  [20]: 


lim 

k-»0 


KQcO 

LK(k) 


lln 

K 


(41 

\k/ 


(18) 


The  substitution  of  (9)  into  (18),  followed  by  the  application  of  (8)  leads  to  the  limiting  behavior  of 
the  normalized  capacitance,  gjvoi  as 

^-»8Dln(2)  +  4Dln(  1 +i)  ,  H^O.  (19) 

Grover’s  expression  (6)  was  not  considered  for  vanishingly  small  H  because  it  is  only  a  function 
of  Din  this  case. 
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ii)  Tube  far  from  the  groundplane 

In  this  region,  the  capacitance  of  the  tube  is  independent  of  the  height  above  the  groundplane 
and 

^ -»a(D)  ,  (20) 

Cq* 

In  (20),  a(D)  is  a  function  determined  from  the  method  of  moments  data. 


Expressions  (19)  and  (20)  can  be  combined  to  yield  a  function  possessing  the  limiting 
behaviors  described  above,  as 


il 


a(D)  +  4D  In 


1  + 


P(D) 


H 


(21) 


In  (21),  the  function  a(D)  replaces  8D  ln(2)  and  P(D)  is  introduced  in  order  to  allow  the  second 
term  to  remain  valid  for  large  values  of  D  and  H.  The  factor  4D  in  the  second  term  has  been 
retained  since  it  is  valid  for  small  values  of  H. 

The  normalized  tubular  capacitance,  C/(eo/),  was  computed  using  the  method  of  moments  for 
the  parameter  range  of  -2.5  <  logio(D)  ^  1  and  -4  <  logio(H)  <  1.  A  plot  of  C/(eo/)  as  a  function 
of  D  and  H  resulting  from  the  computation  in  this  range  is  shown  in  Fig.  9.  From  the  data,  the 
functions  a(D)  and  p(D)  ware  determined  by  nonlinear  regression,  with  the  following  forms: 


P(D)  = 


1  3QD  H-  124D2 
70D(D+2) 


(22) 


The  substitution  of  the  expressions  in  (22)  into  (21)  results  in  an  expression  for  the  normalized 
capacitance  valid  in  the  desired  range,  given  by 


-C.  = 

Eq/ 


+  4Dln 


1  + 


1  3QD  12402]  1 
70HD(D  +  2)  J|  ■ 


(23) 


Equation  (23)  will  be  referred  to  as  the  approximate  capacitance  formula  (ACF).  An  error 
contour  plot  of  (23)  is  shown  in  Fig.  10.  The  plot  indicates  that  the  error  in  the  ACF  is  within 
10%  of  the  method  of  moments  results  over  virtually  the  entire  region  shown.  In  particular,  the 
error  in  the  ACF  is  less  than  3%  for  H  <  0. 1  and  generally  increases  with  H.  The  ACF  is  a  simple 
formula  that  provides  an  accurate  estimate  of  capacitance  over  the  designated  region  of  interest. 
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6.  DISCUSSION  OF  RESULTS 

Table  I  summarizes  the  ranges  where  the  capacitance  formulas  discussed  in  this  report  are 
within  10%  of  the  method  of  moments  results.  The  four  formulas  presented  here  are  seen  to  be 
accurate  for  wide  ranges  in  the  parameters  D  and  H.  In  the  case  of  the  CMA  formula,  the  lower 
boundary  defining  the  region  in  which  the  accuracy  is  better  than  10%  is  given  as  an  approximate 
function  of  H. 

Although  in  practice  monopoles  are  mounted  close  to  the  ground  plane  (H  «  1),  we  observe 
that  the  capacitance  formulas  are  accurate  over  a  wide  range  of  H.  Undoubtedly,  the  use  of  a 
tubular  monopole  far  removed  from  the  ground  is  unlikely  in  most  applications.  In  ordo"  to  gauge 
the  utility  of  these  formulas  for  large  separations  from  the  ground,  it  would  be  instructive  to  do  a 
simple  comparison  against  known  capacitance  formulas  for  tubes  of  arbitrary  diameters  and  lengths 
in  free  space. 

For  the  investigation  of  tubular  monopole  capacitance  formulas  for  large  H,  we  define  a  ratio  R 
as 


R  =  ^=1.6 

Cfs 


(24) 


whCTe  Cg,  and  Cf,  refer  to  the  capacitances  of  a  tube  above  a  ground  plane  and  in  free  space, 
respectively,  and  5  is  the  difference  in  the  capacitances  relative  to  the  free  space  value.  As  H 
becomes  large,  R  -4  1  (i.e.,  8  — >  0).  A  "free  space"  boundary  can  be  defined  as  the  normalized 
height  H  that  corresponds  to  a  pred^ermined  small  value  of  8  (i.e.,  8  «  1).  A  mathematical  fit 
describing  the  boundary  for  the  height-to-length  ratio  Hfsb  along  which  8  =  0.01  is  given  by 


(25) 


applicable  for  D  >  10  .  Expression  (25),  obtained  from  an  approximate  fit  to  the  method  of 
moments  data,  essentially  separates  two  regions;  for  H  <  the  monopole  formulas  are 
appropriate,  while  for  H  >  the  free  space  expressions  apply. 

Expressions  for  the  free  space  capacitance  of  thin  and  thick  tubes  have  been  derived  by  Howe 
[12]  and  Butler  [21],  respectively.  The  expressions  are  given  as  follows: 

CHowe  =  -4?^  ,  D«1  (26) 


/n  _  271  ^£od 

®“““"ln(16D) 


D»  1. 


(27) 
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Butler  [21]  observed  that  the  parameter  range  for  each  expression  can  be  relaxed  without  serious 
error.  He  noted  that  at  D  =  0.25  both  expressions  produce  results  within  4%  of  the  precise  value 
obtained  by  the  method  of  moments,  thereby  allowing  one  the  means  of  computing  the  capacitance 
by  using  the  appropriate  expression  above  or  below  D  =  0.25. 

Figure  1 1  shows  the  regions  of  acceptable  accuracy  of  the  tubular  monopole  capacitance 
formulas  and  the  free  space  tubular  capacitance  formulas.  Some  of  the  fine  detail  has  been 
removed  for  clarity.  The  extended  Grover  formula  (17)  is  omitted  since  its  useful  region  is  close  to 
that  of  Grover.  Note  that  the  Grover  formula  overlaps  with  that  of  Howe  while  the  CMA  result 
approaches  the  approximate  free  space  boundary  in  an  asymptotic  fashion. 

7.  APPLICATION  TO  A  TUBE  OF  NONCIRCULAR  CROSS  SECTION 

The  concept  of  modeling  an  antenna  of  noncircular  cross  section  with  one  of  circular  cross 
section  having  an  effective  radius  that  yields  the  same  input  impedance  and  radiation  pattern  was 
first  introduced  by  Hallen  [22].  This  concept  applies  to  electrically  small  antennas.  If  the 
capacitance  of  the  feed  region  is  negligible,  the  effective  radius  for  a  monopole  of  noncircular  cross 
section  can  be  computed  with  the  aid  of  expressions  from  several  sources  [23-25]  and  used  in  the 
capacitance  expressions  given  in  this  report.  However,  it  must  be  emphasized  that  these  results  are 
approximate. 


8.  SUMMARY 

Formulas  have  been  presented  for  the  computation  of  the  input  capacitance  of  electrically  small 
tubular  monopoles  for  a  wide  range  of  diameters  and  heights  above  ground.  The  regions  of 
validity  for  the  existing  formulas  were  determined  through  comparison  with  precise  results 
obtained  from  a  standard  method  of  moments  code.  Several  regions  lying  between  the  thin  and 
thick  tube  domains  were  identified  where  the  existing  formulas  are  not  applicable.  As  a  result,  an 
approximate  expression  was  constructed,  effectively  bridging  the  gap  between  the  thin  and  thick 
tube  domains.  Taken  as  a  whole,  such  an  assembly  of  formulas  allows  for  the  determination  of 
capacitance  for  a  continuous  range  in  terms  of  the  normalized  variables  D  and  H  for  over  six  orders 
of  magnitude,  with  errors  usually  much  less  than  10%.  For  large  ground  plane  separations,  it  was 
shown  that  the  monopole  formulas  approach  the  results  obtained  for  a  tube  in  free  space. 

For  information  on  various  relevant  topics  including  capacitance  formulas  for  several  tubular 
monopole  feeds  and  the  numerical  convergence  of  the  method  of  moments  algorithm  used  in  this 
investigation,  the  reader  is  referred  to  the  appendices. 
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Figure  1.  (a)  Tubular  monopole. 

(b)  Coordinate  system  used  for  tubular  monopole  analysis. 
Note:  t  «  d,  /. 
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Figure  2.  Typical  methods  for  feeding  tubular  monopoles. 

(a)  Wire  feed. 

(b)  Conical  feed. 
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Figure  4.  Error  contours  for  Grover's  formula  (6). 
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Figure  7.  Error  contours  for  the  conformal  mappi 


TR  10,817 


Figure  8.  Error  contours  for  the  extended  Grover  formula  (17). 
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Figure  9.  Surface  of  the  normalized  tubular  monopole  capacitance  used  to  determine  ACF  (23). 


CONFORMAL  MAPPING  APPROXIMATION 
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Figure  1 1  Regions  of  validity  (error  <  10%)  of  tubular  monopole  capacitance  formulas 
Note:H  =  h//,D  =  d//. 
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Table  I.  Ranges  of  Validity  of  Tubular  Capacitance  Formulas 

(Error  <  10%) 


Formula 

Range  of  H 

Range  of  D 

(  H  =  h  /  /) 

(  D  =  d/O 

Grover 

<  0.0004 

^0.008 

(6) 

0.0004  <  H  <  0.04 

<  0.27  H  0-45 

^  0.04 

^0.35 

Extended  Grover 

<  0.0005 

<  0.007 

(17) 

0.0005  <  H  <  0.1 

<  0.33  VH" 

>  0.1 

<1 

Approximate  Capacitance 

10-4  ^  H  ^  10 

0.003  ^  D  <  10 

Formula  (ACF) 

(23) 

Conformal  Mapping 

>  10-4 

>  2 

Approximation  (CMA) 

“ln(l4.JL) 

(8) 

i  H/ 
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APPENDIX  A 

CAPACITANCE  FORMULAS  FOR  TUBULAR 
MONOPOLE  FEEDS  AND  OTHER  EFFECTS 

1.  Introduction 

In  this  report,  the  monopole  was  considered  as  a  charged  vertical  tube  above  a  ground  plane. 
The  tube  was  considered  to  have  an  infinitesimally  thin  wall  and  no  connecting  wires.  In  actual 
practice,  a  monopole  antenna  has  some  finite  wall  thickness,  however  small,  and  is  fed  in  a 
definitive  fashion.  The  total  input  capacitance  of  a  practical  antenna  must  therefore  account  for 
such  effects.  In  this  appendix,  some  capacitance  expressions  for  these  effects,  suitable  for 
practical  antenna  analysis  and  design,  are  summarized. 

2.  Summary  of  capacitance  formulas 
2.1.  Correction  for  Hnite  wall  thickness 


A  correction  to  account  for  the  finite  thickness  of  the  tubular  monopole  is  developed  by  solving 
the  analogous  problem  of  the  excess  (or  incremental)  capacitance-per-unit  length  of  two  oppositely 
charged  coplanar  plates  of  thickness  t  and  separation  s,  as  shown  in  Fig.  A- 1(a).  The  correction 
AC  is  here  defined  as  the  difference  of  the  per-unit-length  capacitance  C(t,  s)  of  the  coplanar  plate 
structure  with  finitely  thick  plates  and  the  per-unit-length  capacitance  C(0,  s)  of  the  same  structure 
with  plates  having  zero  thickness  : 

AC  =  C(t,.s)-C(0,  s)  .  .  (A-l) 

This  result  is  then  rearranged  to  yield  the  desired  correction  for  the  case  of  a  charged  vertical  plate 
whose  edge  is  parallel  and  over  a  ground  plane,  as  shown  in  Fig.  A-l(b). 

The  plates  in  Fig.  A- 1(a)  extend  to  infinity  in  both  directions  of  the  vertical  plane  and  the 
correction  is  computed  in  the  direction  normal  to  the  page.  The  quantities  C(t,  s)  and  C(0,  s)  are 
determined  through  conformal  mapping.  When  combined  with  (A-l),  an  implicit  solution  results 
and  is  given  by  Cohn  [26]  as 


AC  =  —  Eo  In 
n 

where  k  is  related  to  t/s  by 


2E(k)-k'^K(k) 

2VIc 


(A-2a) 
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and 


t  _(l+k^)K(ki-2E(kj 
^  2[2E(k)-k'^K(k)] 


{A-2b) 


k'=fr^  . 


(A-2c) 


In  (A-2),  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  the  first  and  second  kinds, 
respectively  [19].  A  plot  of  AC  as  a  function  of  t/s  can  be  obtained  by  evaluating  (A-2a)  and  (A- 
2b)  individually,  by  varying  the  modulus  k  from  zero  to  one.  A  plot  constructed  in  this  way  is 
shown  in  Fig.  A-2.  In  the  plot,  the  value  of  the  difference  [AC-(eot/s)]/eo  as  a  function  of  t/s  was 
obtained  by  varying  k  from  0.001  to  0.999,  in  increments  of  0.001.  The  second  term  of  the 
diffeaice,  (eot/s),  is  the  capadtance-per-unit-length  of  the  gap  region  shown  in  Fig.  A-1, 
neglecting  fringing  effects. 

An  inspection  of  Fig.  A-2  shows  that  for  t/s  >  1,  [AC  -  (eot/s)]/eo  approaches  a  constant  and  is 
negligible  for  t/s  ^  0.001.  A  good  approximation  for  (A-2)  that  eliminates  the  intermediate 
computations  involving  k  and  is  applicable  for  all  values  of  t/s  >  0.001,  has  the  following  form: 


(A-3) 


Expression  (A-3)  is  accurate  to  three  significant  figures.  Refaiing  again  to  Fig.  A-l(b),  the 
corresponding  thickness  correction  for  the  vertical  plate  above  a  ground  plane  is  obtained  by 
substituting  2h  =  s  into  (3),  followed  by  a  doubling  of  the  right  hand  side  of  the  expression.  The 
correction  for  the  unsymmetrical  case  thus  becomes 


Expression  (A-4)  is  accurate  to  three  significant  figures  for  t/h  >  0.001.  The  expression  is 
gaieral  enough  so  that  it  can  be  applied  to  a  tubular  monopole  of  arbitrary  cross  section  whai 
multiplied  by  the  mean  perimeter  of  the  tube.  As  an  example,  the  thickness  correction  for  a  tubular 
monopole  with  a  circular  ctoss  section  is  obtained  by  multiplying  the  per-unit  length  correction  (A- 
4)  by  the  mean  circumference  Pmean  of  the  tube,  i.e.. 


Pmean  =  7t(a  +  b)  , 


where  a  and  b  are  the  inna:  and  outs'  radii  of  the  tube,  respectively.  Upon  substitution  of  Pmean 
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into  (A-4),  the  correction  is  given  as 


The  capacitance  correction  given  by  (A-5)  is  added  to  the  capacitance  computed  for  a  thin  walled 
cylinder.  In  the  application  of  an  appropriate  thin-wall  capacitance  expression  to  a  tube  of  arbitrary 
cross  section  and  wall  thickness,  it  is  suggested  that  the  mean  effective  diameter  be  used. 

The  range  of  validity  of  (A-2)  is  given  by  Cohn  [26]  as; 

and  , 

X  *  ^  ^ 

where  s  is  the  gap  distance  between  the  plates  [see  Fig.  A-1  (a)]  and  r  is  the  distance  from  the 
center  point  of  the  configuration  to  the  nearest  extraneous  surface.  An  example  of  an  extraneous 
surface  is  shown  in  Fig.  A-3  for  the  case  of  a  capacitive  iris  in  a  waveguide.  In  this  figure,  the 
distance  r  is  measured  from  the  center  line  of  the  capacitive  iris  to  the  adjoining  waveguide  walls. 
The  range  of  validity  for  the  corresponding  case  of  the  tubular  monopole  is  found  by  setting  2h  =  s 
and  r  =  h+/  and  substituting  them  into  the  above  inequalities.  The  use  of  r  =  h+/  is  made  because 
an  extraneous  surface  (in  the  form  of  a  top  load)  may  be  placed  at  this  distance  to  increase  the  . 
effective  height  of  the  antenna.  The  ranges  of  validity  for  the  corresponding  tubular  monopole  case 
are  then 

—  <  —  and  —  <  —  . 

S  I  3 

In  the  present  application,  the  monopole  is  electrically  small  so  that  h+/  <  A/S,  thereby  satisfying 

the  first  inequality.  The  second  inequality  is  satisfied  in  most  practical  designs. 

2.2.  Thin  wire  feed  capacitance 


The  geometry  of  a  thin  wire  used  to  feed  a  tubular  monopole  is  shown  in  Fig.  A-4.  As  shown 
in  the  figure,  the  wire  is  a  continuation  of  the  center  conductor  of  a  coaxial  cable  and  is  attached  to 
the  lower  end  of  the  tube.  The  capacitance  of  die  feed  wire  is  given  by  King  [16]  as 


where  h  and  a  are  the  height  and  radius  of  the  feed  wire,  respectively.  Expression  (A-6)  is  valid 
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for  h/a  >  75  and  h  <  X/8. 

2.3.  End  correction  for  a  coaxial  line 

In  conventional  transmission  line  theory,  for  a  given  ctoss  section,  it  is  assumed  that  the 
charge-pCT-unit-length  on  each  conductor  is  equal  in  magnitude  but  opposite  in  sign.  The 
capacitance  per-unit-length  derived  undo"  this  assumption  becomes  a  function  of  only  the  cross 
section.  However,  if  the  same  line  is  connected  to  an  antaina,  the  charge  distribution  in  the  region 
of  transition  undergoes  a  change.  This  has  the  effect  of  altering  the  antenna  admittance  measured  at 
some  distance  from  the  antaina,  compared  with  the  admittance  predicted  using  standard 
transmission  line  expressions. 

Figure  A-5  provides  a  view  of  the  transition  from  a  coaxial  cable  to  a  thin  monopole  antenna. 
The  portion  of  the  transmission  line  along  which  the  charge  distribution  differs  from  that  usually 
assumed,  will  be  referred  to  as  the  transition  region.  The  transition  region  extends  from  w  =  0  to 
w  =  d,  where  w  is  defined  in  Fig.  A-5.  The  measured  apparent  admittance  Ya  of  a  monopole 
antenna  is  given  approximately  by  King  [27]  as 

Ya  =  Yo  +  j  CD  C  T  ,  for  b/a  >  1  , 

where 

Ct=  I  [c(w)-Co]dw  . 

Jo 

In  (A-7),  Yo  is  the  ideal  theoretical  admittance,  Ct  is  the  lumped  equivalent  transition  capacitance, 
that  is  placed  at  w=0  (see  Fig.  A-5)  to  account  for  the  disturbance  in  charge  distribution  in  the 
transition  region  of  the  transmission  line.  The  pCT-unit-length  capacitances  c(w)  and  Co  in  (A-8)  are 
for  the  transition  region  and  coaxial  line,  respectively.  Note  that  Ya  ->  Yo  as  b/a  1.  In  (A-8),  d 
is  the  range  over  whidi  the  differaice  c(w)  -  Co  is  significant,  extending  to  d  =  10b.  It  is  seen  from 
(A-7)  that  the  d^ermination  of  the  lumped  transition  capacitance  Ct,  which  has  a  negative  value,  is 
important  in  order  to  accurately  quantify  the  antaina's  input  admittance. 

King  [27]  theoretically  determined  the  normalized  transition  capacitance,  -  Cj/cJo ,  in  the  limits 
of  large  and  small  outer-inner  conductor  ratios  b/a.  A  plot  of  the  approximate  computation,  along 
with  an  extrapolated  curve  joining  the  two  regions  (daived  by  King),  is  shown  in  Fig;  A-6. 
Included  in  this  figure  are  some  measured  data  obtained  by  Hartig  [28].  From  the  figure,  it  can  be 
seal  that  the  measured  values  display  the  same  general  troid  as  the  theoretical  ones  but  are  larga 
for  b/a  >  2.  Because  of  the  approximate  nature  of  the  theory  used  to  predict  the  lumped  transition 
capacitance,  and  the  disaepancy  in  the  scaling  of  Hartig's  measurement  values.  King  suggests 


(A-7) 

(A-8) 
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that  either  curve  may  be  used  to  determine  -  Cx/cob. 

In  the  absence  of  theordical  or  measured  data,  a  useful  approximation  for  Cj  that  yields  values 
intermediate  between  those  of  King  and  Hartig  has  the  form 


Ct  5  -  27teb 


ln(b/a) 

.3+ln^(b/a) 


2  <  b/a  <  30 . 


(A-9) 


Table  A-1  is  a  chart  comparing  (A-9)  against  both  theoretical  and  measured  values. 


2.4.  Capacitance  of  a  conical  feed 


Aside  from  the  thin  wire  feed,  another  common  way  of  feeding  a  monopole  uses  a  cone,  as 
shown  in  Fig.  A-7.  The  classical  method  of  determining  the  capacitance-per-unit  length  of  an 
irfinite  cone  above  a  ground  plane  has  been  derived  by  Schelkunoff  [29]  as: 


2jteo 


(A-10) 


The  complicated  nature  of  the  boundary  conditions  on  the  edges  of  the  truncated  region  prohibits 
an  analytical  treatment  of  the  problem;  recourse  to  numerical  techniques  becomes  a  practical 
alternative.  Toward  this  end,  Wilton  [17]  determined  the  capacitance  of  finite-length  cones  with 
and  without  a  topcap  above  a  ground  plane  using  the  method  of  moments  for  2.5*  <  0o  ^  87,5*. 

The  capacitance  data  from  Wilton's  study  was  used  to  obtain  an  approximate  expression  for  the 
capacitance  of  a  conical  feed.  In  deriving  the  approximate  expression,  an  approach  was  taken 
whereby  the  normalized  capacitances  of  the  infinite  cone  (C/Eo)  and  truncated  cone  (CVeoL)  were 
plotted  against  the  infinite  cone  half  angle  0o,  to  observe  the  deviation  b^ween  the  two  cases. 
While  some  deviation  in  the  capacitance  curves  for  the  two  cases  was  noted,  the  curves  were 
observed  to  be  very  similar  in  shape,  suggesting  that  the  functional  dep«idaice  given  by  (A-10) 
was  applicable.  The  deviation  in  capacitance  between  the  two  cases  was  consequently  corrected  by 
the  substitution  of  an  effective  infinite  cone  half-angle  0o®^^in  place  of  the  original  cone  half-angle 
0o  in  (A-10).  The  approximate  expression  for  the  capacitance  of  a  truncated  cone  was  thus  derived 
with  the  following  form: 


where 


In 


27t£oL 


(A-ll) 
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and 


=  00 


1  + 


90'-  00 


(deg.) 


K  =  74  (cone  with  a  topcap) 

=  76  (cone  without  a  topcap). 


The  capacitances  determined  from  (A-1 1)  are  within  6%  of  the  m^hod  of  moments  data  for  2.5’  ^ 
00  <  87.5’. 

2.5.  Capacitance  of  an  arbitrarily  shaped  plate  above  a  groundplane 

When  the  base  of  a  short  monopole  has  a  plate  attached  to  it,  the  plate  capacitance  along  with 
the  capacitance  of  the  thin  wire  feed  attached  to  the  monopole  must  be  known.  A  plate  of  arbitrary 
shape  with  surface  area  S  and  perimeter  P  situated  above  a  ground  plane,  with  an  interposed 
dielectric  of  relative  pamittivity  e^,  is  shown  in  Fig.  A-8.  Kuester  [18, 30]  has  detornined  the 
capacitance  of  an  arbitrarily  shaped  plate  for  the  cases  where  either  the  narrowest  dimension  Wq  is 
largCT  than  the  substrate  height  h  or  where  the  widest  dimension  Ww  is  smaller  than  the  substrate 
height.  The  expressions  (rearranged  here  for  simplicity)  are  accurate  within  their  regions  of 
validity  and  are  summarized  below. 

1.  Narrowest  plate  dimension  larger  than  the  substrate  height  (Wq  /  h  >  0.5): 


C5^  +  ^[l+ln(2E)  +  Ai(P)  +  e,A2(er)]  ,  (A-12) 


whae 


r 


J 


an  -a^n  (rt/P) 


Po-P 


sin 


JC  (/-/') 


did/'  (A-13a) 


and 


A2  (fir)  =  j-  In  1  -  0.6735  +  0.0788  j^j  ^ 


+  ln(2jc)  .  {A-13b) 


The  variables  used  in  the  line  integral  Ai(P)  are  illustrated  in  Fig.  A-9  for  an  arbitrarily  shaped 
plate.  In  addition.  Table  A-2  lists  values  of  Ai(P)  for  some  common  geometries. 
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2.  Widest  plate  dimension  smaller  than  the  substrate  height  (W^/  h  <  0.5): 


where 


C  = 


27t£o(£r 


arctan 


'hB(E,)- 
.  a*  . 


(A-14) 


(A-15) 


In  (A-14),  ae  is  the  effective  radius  of  the  plate.  Various  methods  has  been  devised  to  compute  Zq, 
with  most  of  them  numerical  [10, 3 1-36].  Among  the  analytical  methods,  those  of  P61ya  and 
Szego  [37]  and  Fabrikant  [38]  are  simple.  The  method  presented  in  [37]  solves  for  the  limits 
within  which  ae  exists,  using  lower  and  upper  bound  radii  (ai  and  au).  In  contrast,  [38]  utilizes  an 
integral  to  estimate  the  effective  radius,  using  the  plate's  centroid  as  the  origin.  The  former  method 
is  easier  to  use  because  the  area  S  and  perimeter  P  of  the  plate  are  the  only  quantities  required.  In 
obtaining  a  close  estimate  of  ae  using  upper  and  lower  bound  radii,  Kuester  [18]  has  suggested 
taking  the  arithmetic-geometric  mean  (AGM)  of  the  two  quantities,  given  by 


ae,AGM  = - 

2K 


(A-16) 


whCTe 


and 


au  =  4- 


(A-17) 

(A-I8) 


In  (A-16),  K[-]  is  the  complete  elliptic  integral  of  the  first  kind  [19].  Effective  radii  for  various 
plate  shapes,  taken  fix)m  [18],  are  givai  in  Table  A-3.  For  plate  shapes  other  than  those  listed  in 
Table  A-3,  expressions  (A- 16)  to  (A-18)  may  be  used  to  estimate  a^. 
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T 


_L 

T 


(b) 


Figure  A-1.  Semi-infinite  coplanar  plate  cross  section. 

(a)  Symmetrical  arrangement; 

(b)  unsymmetrical  equivalent. 
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Figure  A-3.  Lengthwise  view  of  a  capacitive  iris 
waveguide.  The  nearest  extraneous 
surface  in  the  illustration  is  the 
waveguide  wall.  From  Cohn  [26]. 
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Ground  plane 


Coaxial  line 


.  Thin  feed  wire  /  monopole 


2a 


w=0  (dashed  line) 


T 


2b 


w 


Figure  A-5.  Cross  section  of  coaxial  line  /  monopole  transition 

(w=0).  The  quantities  2a  and  2b  denote  the  inner  and 
outer  diameters,  respectively. 
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Surface  area  (S) 


Figure  A-8.  Plate  of  arbitrary  shape  over  a  ground  plane  with 
interposed  dielectric.  Thickness  (or  height)  of 
dielectric  is  h. 
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TABLE  A-l 


Comparison  of  Coaxial  End-Effect  Approximation 


Ct=  -iTteb 


ln{bla) 

3+ln^(bla) 


with  Theory  and  Experiment 


Value  of  -Ct  ,  Picofarads  (pF)* 


b/a 

b  (inch) 

Theory  [27] 

Approx. 

Measured 

[28] 

2.21 

0.156 

0.048 

0.050 

0.048 

4.00 

0.375 

0.127 

0.130 

0.159 

5.32 

0.375 

0.111 

0.115 

0.118 

5.33 

0.500 

0.147 

0.154 

0.223 

7.09 

0.500 

0.123 

0.132 

0.184 

8.00 

0.750 

0.186 

0.184 

0.270 

10.64 

0.750 

0.139 

0.154 

0.223 

18.88 

1.77 

0.266 

0.259 

0.366 

25.11 

1.77 

0.193 

0.221 

0.264 

*  With  Er  =  1. 
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TABLE  A-2 

Values  of  the  Coefficient  Ai(P)  used  with  Capacitance  Approximation 

C  =  1  +  In  (2^]  +  A  i(P)  +  SrA2(er) 

h  7U  I  ^7Zh>  J 

where 

Aiier) s  In  1  -  0.6735  j^j  +  0.0788 +  In  (2k) 

for  Various  Plate  Geometries 

Plate  shape  Perimeter  (P)  Area  (S)  Ai(P) 
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TABLE  A-3  (Cont’d) 


Plate  shape 


Effective  radius, 


0.3961a 


0.2940a 


0.9149a 
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APPENDIX  B 

APPROXIMATION  FOR  THE  RATIO  OF  THE 
COMPLETE  ELLIPTIC  INTEGRALS  K(k')/K(k) 


1.  Background 


The  design  or  analysis  of  planar  transmission  lines,  such  as  striplines  or  coplanar  waveguides, 
frequently  involves  the  computation  of  the  ratio  of  complete  elliptic  integrals  of  the  first  kind, 
K(kO/K(k)  (often  written  simply  as  K'/K),  where 


de 


(B-l) 


The  argument  k  (the  modulus),  a  function  of  the  line  parameters,  numerically  assumes  a  value 
between  zero  and  one.  The  complementary  modulus,  k',  is  related  to  k  by 

k'  =  V 1  -  k2  .  {B-2). 

Expressions  for  the  characteristic  impedance  (or  alternatively,  the  per-unit-length  capacitance  or 
inductance)  of  planar  transmission  lines  are  derived  by  solving  Laplace's  equation  in  two 
dimensions  using  conformal  mapping  [39]. 

A  plot  of  the  ratio  K(kO/K(k)  is  shown  in  Fig.  B-l  as  a  function  of  k  (solid  line).  Also 
included  in  the  figure  are  two  curves  (dotted  lines)  of  the  ratio  in  the  limit  of  small  (k  0)  and 
large  (k  1)  moduli,  given  as  [20] 


and 


lim 

k^O 


'K(k-)' 

L  K(k)  J 


->^ln 

K 


lim 

k-k  1 


Klk-) 

L  K(k)J 


2b|^ 


k' 


(B-3a) 


(B-3b) 


The  expression  in  (B-3b)  is  derived  through  exploitation  of  the  ratio's  antisymmetric  behavior 
about  the  inflection  point,  at  k  =  1/V2.  This  amounts  to  the  substitution  of  k'  in  place  of  k  in  the 
first  expression  (not  proven  here),  followed  by  the  use  of 
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K(k-)  I 

K(k)/K(k') 


(B-3c) 


Techniques  for  evaluating  the  ratio  K(k')/K(k)  are  numerous  and,  depending  on  the 
application,  may  range  from  the  classical  methods  employing  theta  functions  [40],  to  numerical  or 
iterative  schemes  [19, 41-42].  All  of  these  methods  are  rather  efficient  when  properly  implemented 
on  a  computer.  However,  most  of  them  use  several  expressions  or  repetitions  thereof  to  cover  the 
entire  domain  (0  <  k  <  1),  making  them  unattractive  for  manual  computation.  The  need  for  a 
simple  and  accurate  expression,  valid  over  the  desired  domain,  was  recognized.  Moreover,  the 
simplification  of  complicated  analytical  expressions  in  which  the  ratio  KfkO/KOc)  appears  may  be 
made  through  substitution  of  an  accurate  approximation.  These  two  objectives  provided  the 
motivation  to  undertake  this  study. 


2.  Derivation  of  K(k')/K(k)  approximation 


2.1.  Approach 

To  be  useful,  an  approximation  for  K(kO/K(k)  must  yield  accurate  numerical  values  at  three 
critical  points  along  the  domain: 

Value  of  k  Value  of  K(k')/K(k) 

0 

1/VI 
1 

A  first  approximation  for  K(10/K(k)  is  derived  by  equating  the  approximate  and  exact  formulas  for 
the  characteristic  impedance  of  a  transmission  line.  This  method  is  similar  to  the  one  used  by 
Hilberg  [41-43].  An  improvement  in  the  accuracy  of  the  first  approximation  is  then  made  by 
adding  another  term  to  the  function's  argument,  derived  through  the  aid  of  a  nonlinear  regression 
scheme. 

In  contrast,  Hilberg’s  method  of  improvement  utilizes  an  iterative  stereographic  projection 
method  to  obtain  a  sequence  of  approximations  for  KCkO/Kfk).  This  approach  offers  a  very  high 
degree  of  refinement,  though  not  over  the  entire  domain  of  interest.  The  approximations  derived 
from  Hilberg's  method  are  suitable  only  in  the  domain  spanning  fi'om  the  inflection  point  (k  =  1/VI 
)  to  either  end  point  (k  =  0  or  1).  The  solutions  for  the  remaining  portion  of  the  domain  are  found 
by  interchanging  k'  and  k,  followed  by  the  application  of  (B-3c). 


1 

0 
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2.2.  First  approximation 


Consider  a  transmission  line  whose  cross  section  comprises  two  thin  opposing  curved  plates 
along  a  circle  of  diameter  D  and  separated  by  a  gap  s,  as  shown  in  Fig.  B-2.  An  approximate 
expression  for  the  characteristic  impedance  of  this  line  is  given  as  [44] 


Z„  =  2L. 


arccosh 


cot 


(1)1 


,  (ohms) 


(B-4) 


where 


Tjo  =  V|io/eo  =  wave  impedance  of  free  space  ( =  376.7  ohms ), 


Po  =  47t  X  10 '  ^  H/m  (permeability  of  free  space), 

in 

£o  =  8.854  X  10 '  F/m  (permittivity  of  free  space). 

The  exact  expression  for  the  characteristic  impedance  of  the  same  configuration  is  given  as  [43] 

(B-5) 


Zo  =  ^  •  ,  (ohms) 

2  K(k') 


whCTe  the  modulus  k  =  s/D  =  sin  0.  By  equating  (B-4)  and  (B-5)  and  noting  that 


sin  9 


one  obtains  a  first  approximation  for  K(k')/K(k),  of  the  form 


kM  . 


K(k)  ■  Jt 


=  ^  arccosh 


1  +k' 


{B-6) 


The  accuracy  of  (B-6)  was  determined  by  comparing  the  value  given  by  the  right  hand  side  with  a 
precise  value  of  K(kO/K(k),  computed  with  an  accurate  ( relative  error  less  than  2  x  10  *  * ) 
polynomial  approximation  for  K(k)  [19].  The  comparison  revealed  that  (B-6)  is  accurate  to  within 
three  significant  figures  for  k  <  0.2.  HowevCT,  for  0.2  <  k  <  1,  (B-6)  yields  values  that  are 
smaller  than  those  plotted  in  Fig.  B-1,  suggesting  that  the  required  argumait  must  be  larger  than 
( 1  +k')/k  by  some  amount,  to  yield  better  accuracy. 
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2.3.  Improvement  of  the  first  approximation 


The  domain  over  which  (B-6)  is  valid  can  be  considerably  widened  by  adding  a  correction  term 
in  the  argument,  i.e.,  we  let 


K(k) 


=  ■2-  arccosh 
7t 


1  +k' 
k 


+  5 


Figure  B-3  shows  the  value  of  5  (plotted  against  k',  for  clarity)  required  to  yield  the  correct  value 
of  K(k')/K(k)  over  the  entire  domain  of  0  ^  k  ^  1.  The  values  of  5  were  computed  by  equating  the 
left  and  right  hand  sides  of  the  expression  above,  with  K(k')/K(k)  obtained  by  numerical 
integration.  The  simplest  approximation  for  5,  found  by  applying  nonlinear  regression,  is  given  by 


/  1-k^  _  k-v^ 

V  1+k'  4(l+k')  * 


The  substitution  of  the  above  expression  for  5  into  the  postulated  form  for  K(k')/K(k)  leads  to  the 
following  final  result 


K(k'). 


K(k)  ~K 


=  ^  arccosh 


1  +k^ 

^  ^4(1 +k') 


(B-7) 


The  table  below  compares  (B-7)  with  the  exact  values  at  the  critical  points  given  earlier: 

Value  of  K(k')/K(k) 


Value  of  k 

Exact 

From  (B-7) 

0 

00 

00 

1 

0.999  998 

1 

0 

0 

The  relative  error  of  (B-7)  is  shown  in  Fig.  B-4.  The  curve  was  obtained  through  comparison  of 
the  value  given  by  (B-7)  with  the  exact  value  of  K(k0/K(k),  obtained  by  ntimerical  integration  for 
0.15  <  k  ^  1.  The  numerical  integration  results  are  accurate  to  within  10'^^. 
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3.  Discussion  of  results 

Figure  B-4  indicates  that  the  relative  error  of  (B-7)  is  less  than  2  x  10 '  for  0  <  k  <  0.99. 

For  0.99  <  k  <  0.999,  the  relative  error  increases  in  a  somewhat  oscillatory  fashion,  reaching 
3.5  X  10  ■  ^  at  k  =  0.999.  For  k  >  0.999,  the  relative  error  increases  further  but  quickly  drops  to 
zero  when  k  =  1.  An  understanding  of  the  behavior  of  (B-7)  in  the  region  0.999  <  k  <  1  can  be 
facilitated  with  the  observations  that  follow. 

In  the  limit  of  large  moduli,  (B-7)  reduces  to  the  following  form: 

I  7 

lim  (  —  arccosh 
k^i  I 

In  contrast,  from  (B-3b), 

lim 

1 

U'l 

Both  (B-8a)  and  (B-8b)  numerically  approach  zero  as  k  — »  1  and  are  slowly  varying  in  that  region. 
In  order  to  clearly  distinguish  the  error  incurred  by  the  use  of  (B-7),  it  is  useful  to  compute  the 
ratio  of  the  values  of  (B-8a)  and  (B-8b).  In  doing  so,  one  obtains  the  following  result: 

A  plot  of  (B-9)  is  shown  in  Fig.  B-5.  From  the  figure,  it  can  be  seen  that  (B-7)  is  slightly  larger 
than  (B-3b)  at  logio(k')  =  -2.87  ( k'  =  1.34  x  10  ‘  ^ )  and  is  smaller  elsewhere.  Table  B-1  shows 
how  (B-7)  compares  with  the  exact  values  of  K(kO/K(k),  for  moduli  very  close  to  unity. 

As  shown  in  Table  B-1,  approximation  (B-7)  is  smaller  in  value  than  the  exact  result, 
accounting  for  the  rapid  increase  in  error.  In  practical  work,  howeva*,  values  of  k  such  as  shown 
in  the  table  are  uncommon.  Although  it  can  be  argued  that  (B-3b)  can  be  used  to  rigorously  extend 
the  allowable  range  of  k,  it  is  not  necessary.  In  general,  transmission  line  parameters  derived  from 
such  near-unity  moduli  may  be  eitha-  too  small  or  too  large  to  be  useful  in  practical  applications. 
Therefore,  (B-7)  is  an  expression  that  may  be  used  for  most  practical  computations. 


kM 

LK(k)J 


7C 


(B-8b) 


1  +k' 


k7^ 


4(l+k')J)  ^ 


(B-8a) 
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Figure  B-1.  Ratio  of  the  complete  elliptic  integrals  of  the  first  kind,  K(k')/K(k);  solid  line. 

Shown  in  dotted  lines  are  the  limiting  values  of  KfkO/KCk),  obtained  from  (B-3). 
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TABLE  B-l 


Comparison  of  the  Approximation 


K(k) 


=  2.  arccosh 
n 


1  +  k'  ^  kW 
k  4  (  1  +  A:' ) 


with  Exact  Values,  for  Modulus  k  Close  to  Unity 


Value  of  K(k')/K(k) 


Value  of  k' 

Value  of  k 

Exact 

Approx. 

f 

O 

0.999  950 

0.262  172 

0.265  632 

10-“ 

0.999  999  995 

0.148  235 

0.142  333 

10'* 

l-(5xl0*^^) 

0.103  330 

0.080  003 

o 

1 

oo 

l-(5xl0‘^'^) 

0.079  305 

0.045  007 

10- 10 

l-(5xl0'^^) 

0.064  345 

0.025  312 

0 

1 

»  •  * 

0 

0 
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APPENDIX  C 

NUMERICAL  CONVERGENCE  OF  THE 
METHOD  OF  MOMENTS  ALGORITHM 

1.  Introduction 

The  method  of  moments  code  that  was  used  in  this  investigation  is  based  on  the  solution  of  the 
potaitial  integral  equation  (3)  for  the  unknown  charge  density  q(z).  The  details  of  this  program  are 
provided  in  [15],  Because  the  method  of  momaits  data  are  used  as  a  basis  for  evaluation  of  the 
capacitance  formulas  presorted  in  this  report,  it  is  important  that  a  sufficient  number  of  basis 
functions  are  used  in  ordo  to  insure  convergence  of  the  tubular  monopole  capacitance  to  a  given 
number  of  significant  figures.  This  concern  prompted  a  detailed  numerical  study  of  the 
condition(s)  required  to  insure  convergence.  This  appendix  provides  a  summary  of  this  study. 

2.  Results 

Convergence  tests  were  performed  to  determine  the  required  number  of  basis  functions  for 
diffoent  ranges  of  logio(D)  and  logio(H).  Figure  C-1  shows  the  normalized  capacitance  as  a 
function  of  the  number  of  basis  functions  for  the  case  logjoCD)  =  0  and  log^oCH)  =  -1-  The  graph 
indicates  that  the  data  converges  to  three  significant  figures  with  less  than  1500  basis  functions. 
Figure  C-2  illustrates  a  similar  plot  for  the  same  thickness  tube  with  logio(H)  =  -3.  In  this  case, 
approximately  2500  basis  functions  are  required  before  convCTgence  to  three  significant  figures  is 
attained.  In  general,  the  convCTgaice  of  the  tubular  monopole  capacitance  was  found  to  be  slower 
both  for  thicker  tubes  and  for  smaller  separations  above  the  ground  plane. 

Figure  C-3  shows  the  number  of  basis  functions  used  in  the  method  of  moments  algorithm  for 
diffCTent  ranges  of  tube  thicknesses  and  separations  above  the  ground  plane.  The  basis  functions 
listed  in  Fig.  C-3  provide  convergence  to  at  least  three  significant  figures.  Because  of  the 
variability  in  convergence  with  tube  thickness.  Fig.  C-3  is  divided  into  three  regions,  i.e.,  thin, 
medium,  and  fat.  Since  the  convagence  of  the  tubular  monopole  capacitance  slows  considerably 
for  logio(H)  <  -3,  the  numbo’  of  basis  functions  N  used  in  this  report  is  given  as  a  linear  function 
of  logio(H)  for  each  tube-thickness  region.  In  the  thin-tube  region  (-3  <  logjoCD)  ^  -2),  we  have 

N  =  1900  -  6200  [  logio(H)  +  3.5  ]  ,  -4  <  logio(H)  <  -3.5  .  (C-1) 

Similarly,  in  the  medium-tube  region  (-2  <  logio(D)  <  0),  the  number  of  basis  functions  is  given 
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by 


N=  1900-4100[logio(H)  +  3]  ,  -4 <  log  10(H) < -3  , 


while  in  the  fat-tube  region  (0  <  log  10(D)  <  3),  we  have 


N  =  2500  -  5500  [  logio(H)  +  3  ]  ,  -4  <  log  10(H)  <  -3  •. 


(C-2) 


(C-3| 


Note  that  at  the  smallest  ground-plane  separation  (i.e.,  log  10(H)  =  -4),  5000  basis  functions  are 
required  for  a  thin  tube  in  contrast  to  8000  basis  functions  for  a  fat  tube. 

As  an  additional  check  of  the  algorithm,  the  moment-method  capacitance  data  was  examined  for 
tubular  monopoles  of  different  radii  as  a  function  of  the  height  above  the  ground  plane.  Figures  C- 
4  and  C-5  provide  plots  of  the  normalized  capacitance  as  a  function  of  logio(H)  for  logio(D)  =  -1 
and  logio(D)  =  3,  respectively.  In  both  cases,  the  curves  show  that  the  capacitance  increases 
almost  linearly  with  decreasing  log  10(H).  Monopoles  of  different  diameters  were  observed  to 
exhibit  a  similar  behavior. 
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Normalized  capacitance  of  a  tubular  monopole  versus  the  number  of  basis  functions 
for  log  (D)  =  0  and  log  (H)  =  -1  as  determined  by  the  m^hod  of  moments  algorithm 
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Figure  C-2.  Normalized  capacitance  of  a  tubular  monopole  versus  the  number  of  basis  functions 

for  log  (D)  =  0  and  log  ^(H)  =  -3  as  determined  by  the  method  of  moments  algorithm. 


Refer  to  Eqs.  (C-1),  (C-2),  and  (C-3) 
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log(H) 

Figure  C-3.  Number  of  basis  functions  required  to  insure  a  minimum  of  three  significant  figure  accuracy 
in  the  method  of  moments  algorithm  used  for  computing  tubular  monopole  capacitance. 
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Figure  C-5.  Variation  of  the  normalized  tubular  monopole  capacitance 
with  height  for  logj^(D)  =  3. 
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